Controlling spatiotemporal chaos in excitable media using an array of control points 
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The dynamics of activation waves in excitable media can give rise to spiral turbulence, the resulting 
spatiotemporal chaos being associated with empirical biological phenomena such as life-threatening 
disturbances in the natural rhythm of the heart. In this paper, we propose a spatially extended but 
non-global scheme using an array of control points for terminating such spatiotemporally chaotic 
excitations. A low-amplitude control signal is applied sequentially at each point on the array, 
resulting in a traveling wave of excitation in the underlying medium which drives away the turbulent 
activity. Our method is robust even in the presence of significant heterogeneities in the medium, 
which have often been an impediment to the success of other control schemes. 

PACS numbers: 05.45.Gg,87.18.Hf,87.19.Hh 



Excitable media are a class of models for a wide range 
of physical, chemical and biological systems that show 
spontaneous formation of spatial patterns, such as spiral 
waves [1]. These spiral waves, under certain conditions, 
may become unstable and break up, giving rise to spa- 
tiotemporal chaos. Such phenomena have been seen to 
occur in a large variety of natural systems described by 
excitable media, e.g., through interactions between chem- 
ical concentration waves in the Belusov-Zhabotinsky re- 
action [2] and CO oxidation on Pt(llO) surface [3], cyclic 
AMP waves involved in the organization of multicellular 
morphogenesis in Dictyostelium Q , intracellular calcium 
waves in Xenopus oocytes [H[ , intercellular calcium waves 
in brain slices [y, lZ[ , and propagating electrical activity 
in the pregnant uterus Q and cardiac tissue 

0,113. The 

last example is of particular importance as spatiotem- 
porally varying patterns of excitation have been impli- 
cated in clinically significant disturbances of the natural 
rhythm of the heart, which are termed as arrhythmias. 
In fibrillation, a potentially life-threatening arrhythmia, 
there is complete loss of coordination between different 
regions of the heart due to spatially extended chaotic 
activity The resulting cessation in the mechanical 

pumping action necessary for blood circulation, can lead 
to death within minutes unless the condition is termi- 
nated promptly. Conventional defibrillation is done by 
applying a large electrical shock across the heart, which is 
not only extremely painful, but can also damage cardiac 
tissue. Thus, devising a low-amplitude control method 
for spatiotemporal chaos in excitable media is an excit- 
ing challenge, as well as, having potential clinical rele- 
vance EEEl. 



Low-amplitude control of chaos in excitable media 
needs to take into account certain special features of 
such systems, which hinder the straightforward applica- 
tion of methods successful in other chaotic systems. In 
particular, there exists a threshold for the stimulation, 
only on exceeding which excitation occurs as indicated 
by a characteristic action potential Thus, the control 
signal amplitude needs to be above a lower bound to ef- 



fect discernible changes in the medium, ruling out simple 
proportional feedback schemes. The other complicating 
factor is the existence of a recovery period immediately 
following excitation, during which the medium will not 
respond to a normal supra-threshold stimulus. This im- 
plies that two excitation wavefronts will annihilate on 
collision, as neither can penetrate the other, their wave- 
backs being in the recovery phase. This recovery period 
is distinguished into absolute, during which the cell is in- 
sensitive to any perturbation, and relative, when a larger 
than normal stimulus can excite the cell, the amplitude 
depending on the exact phase of recovery. Therefore, the 
amplitude and timing of the control signal needs to be 
appropriately chosen, so that it results in a response from 
the medium. Note that, the excited state is metastable, 
and the cell eventually recovers to the resting state asso- 
ciated in different biological systems with a characteristic 
resting transmembrane potential (p± —84 mV for cardiac 
myocyte cells). Thus, the control of spatiotemporal chaos 
in excitable media is essentially a problem of synchroniz- 
ing the excitation phase of every cell, so that the entire 
system returns to the resting state, resulting in the ter- 
mination of all activity. 

Chaos control schemes for excitable media may be 
broadly classified into [15[: (i) global, where the con- 
trol signal is applied to all points of the system, (ii) lo- 
cal, where only a small localized region of the system is 
subject to control, and (hi) spatially extended but non- 
global schemes. Non-global methods use less power and 
also are relatively easier to implement practically, need- 
ing fewer control points. However, strictly local control 
methods almost always involve very high-frequency stim- 
ulation [16[ , that can by itself lead to reentrant waves in 



the presence of inhomogeneities Moreover, the 

effect of local stimulation at a point can affect the rest 
of the system only through diffusion. As wavefronts an- 
nihilate on collision, control-induced waves are restricted 
to the local neighborhood of the stimulation point dur- 
ing spiral turbulence, with the existing excited fragments 
closer to the control point shielding chaotic activity fur- 



2 



ther away. By using a spatially extended but non-global 
scheme [19j one can potentially avoid these drawbacks. 
In this paper, we terminate chaos in excitable media by 
applying spatiotemporally varying stimulation along an 
array of control points. The control signal appears to 
propagate along the array, triggering an excitation wave- 
front in the underlying medium, that is regenerated after 
each collision with chaotic fragments and in the process 
eliminating all existing activity. Stimulating each point 
once (or at most, twice) is seen to successfully control 
chaos in almost all instances. Although an array of con- 
trol points have been used earlier to prevent the breakup 
of a single spiral [20[ , to the best of our knowledge this is 
the first instance showing control of fully developed spa- 
tiotemporal chaos in excitable media using only a finite 
number of control points without repeated stimulations 
at high frequency. 

The spatiotemporal dynamics of excitation in several 
biological systems can be described by: 
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where V (mV) is the transmembrane potential, C m 
= 1 /iF cm -2 is the transmembrane capacitance, D 
(cm 2 s _1 ) is the diffusion constant, Ii on (/^A cm -2 ) is 
the transmembrane ionic current density and I ext (x, y, t) 
is the space- and time-dependent external stimulus cur- 
rent density applied for control on a 2-dimensional sur- 
face. For the specific functional form of I ion , we used the 



Luo-Rudy I (LR1) action potential model [2l[, where, 
it is assumed to be composed of six distinct currents, 
each of them being determined by several time-dependent 
ion-channel gating variables £ whose time-evolution is 
described by differential equations ^ = The 
parameters in these equations are the steady-state val- 
ues of £, = a^/(a^ + and the time constants, 
Tf = l/(o6£ + /%), which are governed by the voltage- 
dependent rate constants for the opening and closing of 
the channels, and themselves complicated func- 
tions of V. In order to verify the model independence 
of our results and to carry out three-dimensional simu- 
lations, we have also used a simpler description of the 
action potential, as given by Panfilov (PV) [22|, |23| : 
hon = f(V)—g, where / is a piece wise linear approxima- 
tion of a cubic function and g is an effective membrane 
conductance evolving with time as ^| = e(V,g)(kV — g). 
The time-constant e is a function of both V and g and 
the parameters used are same as that in Ref. [19]. The 
models are solved using a forward-Euler scheme, the sys- 
tem being discretized on a spatial grid with spacing Sx 
(=0.0225 cm for LR1, = 0.05 cm for PV). The simula- 
tion domain is a square lattice of L x L points in two 
dimensions or a cuboid with L x L x L z points in three 
dimensions. For the LR1 simulations, L = 400, while for 
PV, L = 256 (for 2-d) and L = 128, L z = 8 (for 3-d). The 
standard five-point and seven-point difference stencils are 



used for the Laplacian in two and three dimensions, re- 
spectively. The time step for integration is chosen to be 
St = 0.01 ms (for LR1) and = 0.11 ms (for PV). No-flux 
boundary conditions are implemented at the edges of the 
simulation domain. The initial spatiotemporally chaotic 
state is obtained by creating a broken wavefront which 
evolves into a spiral wave and is then allowed to become 
unstable, eventually breaking up into multiple wavelets 
(in LR1 this process takes, on average, 200 ms). 

We now focus on the control term I ext (x,y,t). For a 
2-d domain of size L x L we consider 



I ext = I(x, y, t)S(x - md)S(y - nd), 



(2) 



where, the delta function is defined as S(z) = 1 if z = 0, 
and = 0, otherwise, d is the spatial interval between 
points in an array where the control signal is applied and 
m, n are integers in the interval [0, L/d] . The current den- 
sity I(x,y,t) = J for t £ [^/x 2 + y 2 /v, x 2 + y 2 /v) + 
r], and = otherwise, corresponds to a rectangular con- 
trol pulse of amplitude Io of duration r that is travelling 
with velocity v. At the onset of control {t = 0), the point 
at (0, 0) is stimulated, followed a short duration later by 
the points at (0, d) and (d, 0), and this process continues 
as the control pulse proceeds like a traveling wave across 
the array. At each control point, the stimulation may 
excite the underlying region depending on its recovery 
phase. An excited region can in turn spread the effect 
of the stimulation to the surrounding regions through 
diffusion. Fig. [T] shows this process of secondary wave 
generation at the stimulated control points, creating a 
sustained excitation wavefront. Note that, any portion 
of this stimulated wavefront in the medium that is broken 
through collision with chaotic fragments, can be regen- 
erated by subsequent control points. Hence, there is ef- 
fectively an unbroken wavefront that travels through the 
medium, sweeping away the spatiotemporal chaos and 
leaving the system in a recovering state. As each control 
point needs to impose order over a region of size ~ d 2 to 
eliminate spatiotemporal chaos, this highlights the criti- 
cal role of d. Indeed, for large v, the method approaches 
global control as d — > 0, while for d ~ L the control 
stimulation is confined to a single, localized region. This 
implies that as d increases, terminating chaos becomes 
increasingly difficult. For example, in LR1 model, chaos 
control fails for d > 13. 

Before reporting the simulation results, we consider 
the role played by the traveling wave nature of the con- 
trol pulse. The success of the proposed method depends 
on the control signal velocity, v, relative to the velocity 
of the excitation wavefront propagating via diffusion, c 
(= 60.75 cm s _1 for LR-1, = 50.9 cm s _1 for PV, for 
the parameters used here). When v — > 0, the control 
method reduces to a local scheme regardless of d, as the 
effect of the external stimulation can only propagate in 
the system via diffusion. On the other hand, when v is 
very large, all the control points are stimulated almost 



FIG. 1: (left) Schematic diagram and (right) pseudo- gray- 
scale plot of the transmembrane potential V for the two- 
dimensional LR1 model, showing the propagation of the 
control-induced excitation wavefront. The control signal trav- 
eling with velocity v stimulates at time t\ a column of points 
(spaced d apart) from which secondary excitations traveling 
with velocity c are generated. At time £2 the control signal 
stimulates the next column of points, which the induced ex- 
citation wave may or may not have reached depending on the 
relative values of v and c [v/c = 4 in (right)]. 



simultaneously. While the traveling wave nature of the 
control pulse allows propagation of stimulation indepen- 
dent of diffusion through the excitable medium, for v ~ c 
the excitation propagating by diffusion reinforces the ex- 
ternal stimulation at each control point. Further, for a 
control signal propagating with a finite velocity (thus en- 
gaging only a few points at any given time), the energy 
applied per unit time to the medium is much lower than 
that for simultaneous stimulation of all the control points 
(i.e., v — > 00). 

For a system undergoing chaotic activity, the medium 
will at any time be at an extremely heterogeneous state, 
with certain regions excited and other regions partially 
or fully recovered. The vital condition for successful ter- 
mination of chaos is that after the passage of the control- 
stimulated wave there should not remain any unexcited 
region which is partially recovered and can be subse- 
quently activated by diffusion from a decaying excitation 
front. This places a lower bound on the control signal 
parameters, i.e., the signal amplitude io and its duration 
r. If either is decreased below this bound, the external 
stimulation is unable to excite certain partially recovered 
regions. If these regions have neighboring chaotic frag- 
ments, whose activity is slowly decaying after collision 
with the control-induced wave, then, there will be a dif- 
fusion current from the latter. Depending on the phase 
of recovery, this may be sufficient to stimulate activity 
in the partially recovered regions, thereby re-initiating 
spatiotemporal chaos after the control signal has passed 
through. 

To understand in detail the lower bound on the exter- 
nal stimulation parameters io and r, we first look at the 
condition for exciting a completely homogeneous medium 
in the resting state. The stimulation at each point must 
exceed the local threshold in order to generate an action 
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FIG. 2: Strength-duration curves for LR1 model (in a two- 
dimensional domain with L = 400) when control stimula- 
tion is applied on (left) quiescent homogeneous medium and 
(right) medium with existing excitation activity, either a sin- 
gle spiral or spatiotemporal chaos. In both figures, control 
is applied over a grid of points which are spaced apart by 
d — 10. Different curves correspond to different control signal 
velocities v, relative to c, the excitation wavefront velocity in 
the medium. For values of io and r above the curves, the 
external stimulation results in (left) excitation of the domain 
or (right) control of existing activity. 



potential. This could be achieved either directly through 
an external current I ext or indirectly through diffusion 
from a neighboring excited region. Fig. [2] (left) shows 
the result of applying control signals with different Io 
and r at points which are spaced a distance d apart. 
The resulting strength- duration curve 0, HH indicates 
that the response of the system is not sensitively depen- 
dent on the propagation velocity v of the control signal 
along the grid. As d decreases, excitation is possible at 
lower values of Io and r, the minimum being for the case 
when all points are subject to direct external stimulation 
(d — > 0). This is because the entire applied current Io at 
any point is used to raise its state above the threshold, no 
part being lost to neighboring regions through diffusion. 

For systems with existing activity, such as self- 
sustaining spiral waves or spatiotemporal chaos, the re- 
gions in the relative recovery period can be excited by 
stimuli larger than that needed for a fully recovered 
medium. Hence, the strength-duration curve for control 
of such a system will shift towards higher values of Io and 
r (Fig. [21 right). We observe that the minimum external 
stimulus required for control does not vary significantly 
if the medium is undergoing fully developed spatiotem- 
poral chaos as opposed to having a single spiral wave. 
Further, the velocity of the control signal along the array 
is not critical to the success in eliminating existing activ- 
ity, provided v is not significantly smaller than c. Note 
that, if v is sufficiently small or d increases beyond a 
critical value, the control fails, as the effect of the control 
signal is confined to the region immediately surrounding 
the stimulated point. 

Fig. [3] shows the successful termination of spatiotem- 
poral chaos in the LR1 model using a control signal that 
travels across the 2-dimensional domain while exciting 
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T = 90 ms 




FIG. 3: Pseudo- gray- scale plots of the transmembrane poten- 
tial V for the two-dimensional LR1 model (L = 400) showing 
the elimination of all chaotic activity within 300 ms after ini- 
tiation of control. A single wave of control stimulus (Jo = 75 
fiA I cm 2 , r = 3 ms) begins at the top left corner (T = 0) and 
travels across the domain with velocity v — c. This results 
in stimulating the region around the control points (spaced 
apart by d — 10) in a sequential manner, creating a stimu- 
lated wavefront seen as an arc consisting of excited points in 
the panels above. 



points in the medium that are spaced apart by d = 10. 
We verified that the control scheme is not model depen- 
dent by using it to eliminate spatiotemporal chaotic ac- 
tivity in the PV model. As most systems in reality have 
depth, it is crucial to verify that the method is success- 
ful in controlling chaos in a 3-dimensional domain, even 
when the external stimulus is applied only on one sur- 
face. This latter restriction follows from the fact that, in 
most practical situations it may not be possible (or de- 
sirable) to penetrate the medium physically in order to 
apply control signals inside the bulk. We confirmed that 
our method works in thin slices of excitable media of size 
L x L x L z (L z <C L), when the array of control points 
is placed on one of the L x L surfaces (Fig. [4]). Even in 
cases where a single control-stimulated wave across the 
medium is unable to terminate all activity, we notice that 
it results in driving the chaotic activity further towards 
the boundaries and away from the origin of control stim- 
ulation. Thus, using multiple waves through application 
of control signals at intervals which are larger than the 
recovery period of the medium, the chaos in the bulk of 
3-dimensional systems is successfully terminated. 

We have checked that small distortions in the regular 
array of control points does not result in the failure of our 
method. Similarly, starting the control signal at different 
points of origin (and indeed, using a planar wave rather 
than a curved wave) does not affect the efficacy of the 




FIG. 4: The effect of applying control on a single surface of 
a 3-dimensional domain using the PV model (L = 128 and 
L z = 8). Two waves of control signal (Jo = 10, r = 16.5 ms, 
v — c) are applied 165 ms apart resulting in termination of 
spatiotemporal chaos. The interval between control points, 
d — ► 0. The panels show isosurface plots at T = 242 ms (left), 
308 ms (center) and 341 ms (right). 
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FIG. 5: Pseudo- gray-scale plots of the transmembrane poten- 
tial V for the two-dimensional PV model (L = 256) showing 
chaos control in the presence of a conduction inhomogene- 
ity (of size 110 x 110, indicated by the broken lines). Inside 
this region, Di n homogen — 0.0LD, the diffusion constant of the 
rest of the medium. A single pulse of control stimulus (Jo = 8, 
r = 17 ms, v = c) is applied over an array of control points 
spaced apart by d = 8, resulting in termination of all activity 
by T = 350 ms. 



scheme. Further, our method is robust in the presence 
of conduction inhomogeneities (such as inexcitable obsta- 
cles) that tend to destabilise local control schemes. Fig. [5] 
shows successful control of chaos when the medium con- 
tains a large region of slow conduction, i.e., an extremely 
small value of D compared to the rest of the medium. 

In this paper, we have presented a novel control scheme 
involving external stimulation applied over an array of 
points, that is successful in terminating spatiotemporal 
chaos in both simplified as well as realistic models of bi- 
ological excitable media. The control signal amplitude is 
varied both spatially and temporally, such that it appears 
as a propagating wave along the array of control points. 
This results in a stimulated wavefront in the excitable 
medium, that, depending on the propagation velocity of 
the control signal and the space interval between con- 
trol points, eliminates all existing activity. Our method 
requires very low-amplitude control currents applied for 
short durations at a finite number of points, each point 
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being stimulated once (or at most, twice) in most situ- 
ations. Further, it is successful in terminating chaos in 
the bulk of a three-dimensional medium even when ap- 
plied only on one surface. The use of significantly lower 
number of control points than that necessary for global 
control methods, makes the proposed scheme more suit- 
able for practical implementation. 

We thank the IMSc Complex Systems Project and 
IFCPAR (Project 3404-4) for financial support. 
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